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Abstract —  Douglas-KoII-Hess  (DKH)  second-order  relativistic  scalar  approach  was  used  to  investigate 
piezoeletricity  in  gadolinium  ferrite  (GdFeOj).  To  adequately  represent  the  polyatomic  environment 
studied  -  (24sl3p),  (29sl7pl2d)  and  (32s22pl6dl0fl-  basis  sets  were  built  for  the  atoms  O  (3P),  Fe  (5D), 
and  Gd  (9D),  then  contracted  to  [4s2p],  [12s6p5d],  and  [19sl2p8d4f],  respectively.  The  qualifying  of  the 
contracted  basis  sets  for  GdFeOi  crystal  studies  was  conducted  in  three  moments,  namely:  quality 
evaluation  in  molecular  calculations,  made  in  1 FeOI+  and  'GdOl+  molecular  fragments;  the  choice  of 
polarization  function  used  in  the  [4s2p]  basis  set  for  O  (3P)  atom;  the  choice  of  diffuse  functions  used  in  the 
[12s6p5d]  and  [19sl2p8d4f]  basis  sets  for  Fe  (’D)  and  Gd  (9D)  atoms,  respectively.  The  qualified 
contracted  basis  sets  gave  rise  to  the  molecular  [4s2pld]/[13s7p6]/[20sl3p9d5f]  basis  set  which  was  then 
used  to  describe  geometric  parameters  of  the  GdFeOs  crystal.  The  good  performance  of  the 
[4s2pld]/[13s7p6d]/[20sl3p9d5f]  basis  set  in  describing  the  geometry  of  the  material  of  interest  led  to 
calculations  of  the  material  properties:  total  relativistic  energy,  dipole  moment  and  total  atomic  charge. 
The  analysis  of  the  results  for  these  properties  showed  that  the  possible  piezoelectricity  in  GdFeO}  can  be 
caused  by  electrostatic  interactions  of  its  atoms. 
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I.  INTRODUCTION 

Since  the  discovery  of  the  piezoelectric  effect  on 
Barium  Titanate  (BaTiO.i)  ceramics  by  Roberts  in  1947  [1], 
a  large  number  of  perovskite  oxides  presenting  this 
property  have  been  obtained.  Perovskites  are  of  great 
interest  in  materials  science  because  they  are  relatively 
simple  crystalline  structures  and  exhibit  many  properties 
such  as  electrical,  magnetic,  optical,  and  catalytic 
properties  among  others  [2-7]. 

Certain  crystalline  materials  have  the  ability  to 
develop  an  electric  charge  propotional  to  a  mechanical 
stress  called  piezoelectricity.  It  has  been  realized  that 
materials  showing  this  phenomenon  must  also  show  the 
converse,  a  geometric  strain  (deformation)  proportional  to 
an  applied  voltage  [8]. 

www.iiaers.com 


The  perovskite  structure  is  anetwork  of  corner 
linked  oxygen  octahedra,  with  the  smaller  cations  filling 
the  octahedral  holes  and  the  large  cations  filling  the 
dodecahedral  holes  [9],  The  piezoeletric  properties  in 
perovskite  structure  result  from  uncentersymmetric 
characteristics  [10], 

We  have  reported  in  the  literature  computational 
investigations  on  piezolectric  properties  in  perovskites.  In 
these  studies,  the  importance  of  the  set  of  atomic  bases 
developed  exclusively  to  adequately  represent  the 
polyatomic  environment  has  been  pointed  out  [11,  12].  The 
details  of  this  subject  can  be  found  elsewhere  [13], 

Results  from  the  theory  canhelp  experimentalists 
better  design  their  experiments  to  rationalize  the  use  of 
time  and  resources  to  study  a  system  under  investigation. 
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The  aim  of  this  study  is  to  provide  some  insight  into  the 
investigation  of  gadolinium  ferrite  (GdFeCh) 
piezoelectricity  through  the  use  of  quantum  mechanics. 

The  very- well  documented  generator  coordinate 
Hartree-Fock  (GCHF)  method  [13-16]  was  used  to 
build(24sl3p),  (29sl7pl2d),  and  (32s22pl6dl0f)  basis  sets 
for  O  (3P),  Fe  (5D),  and  Gd  (9D)  atoms,  respectively, 
whichwere  contracted  to  [4s2p],  [12s6p5d],  and 

[19sl2p8d4f],  evaluated  in  molecular  calculations,  and 
then  after  supplementation  of  these  sets  with  polarization 
functions  and  diffuse  functions,  to  adequately  represent  the 
polydomic  environment  of  the  GdFeOs  crystal,  through  a 
good  description  of  the  total  energy  and  the  energies  of  the 
F10M0(highest  occupied  molecular  orbital)  and  FlOMO-1 
orbitals  (one  level  below  to  highest  occupied  molecular 
orbital).  Lastly,  contracted  and  supplemented  basis  sets 
with  polarization  function  and  diffuse  functions  were  used 
in  the  quantum-mechanical  study  of  piezoelectricity  in 
GdFe03  at  DKH  second-order  relativistic  scalar  level  [17- 
19]. 

II.  COMPUTATIONAL  METHODOLOGY 

Details  about  the  formalism  as  well  as  the 
strategies  used  by  GCHF  approach  in  the  construction  of 
basis  sets  are  given  in  Refs.  [14-16],  In  this  section,  we 
describe  the  procedures  used  in  constructing  the  contracted 
basis  sets  and  qualifying  the  contracted  basis  sets  for  the 
study  of  piezoelectricity  in  GdFeOj. 

Full  details  about  all  wave  functions  generated  in 
this  work  are  available  by  e-mailing  to:  ciriaco@ufpa.br. 

The  contractions  of  the(24sl3p)/[4s2p], 
(29s  1 7p  1 2d)/[  1 2s6p5d] ,  and  (32s22pl6dl0f)/[19sl2p8d4f] 
basis  sets  were  carried  out  with  the  segmented  contraction 
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scheme  [20]  as  follows: [4s2p]  (16,  4,  3,  1/12,1)); 
[12s6p5d]  (10,  3,  2,  1,  1,  1,  2,  2,  2,  2,  2,  1/5,  5,  3,  2,  1,1/8, 
1,  1,  1,  1);  [19sl2p8d4f]  (11,  2,  2,  2,  1,  1,  1,  1,  1,  1,  1,  1,  1, 
1,  1,  1,  1,  1,  1/9,  3,  1  ,1,  1,  1,  1,  1,  1,  1,  1,  1/6,  4,  1,  1,  1,  1, 
1,  1/7,  1,  1,  1). 

The  qualification  procedure  of  the  contracted  basis 
sets  was  performed  in  three  moments,  as  follows:  quality 
evaluation  of  the  contracted  basis  sets  in  molecular 
calculations  with  the'FeO!+  and  'Gd01+  molecular 
fragments;  the  choice  of  polarization  function  used  in  the 
contracted  basis  set  for  the  0(3P)  atom;  the  choice  of 
diffuse  functions  used  in  the  contracted  basis  sets  for  the 
Fe  (5D)  and  Gd  (9D)  atoms.  These  three  moments  will  be 
detailed  below: 

DKH  second-order  scalarrelativistic  [17-19] 
calculations  were  performed  for  the  TRE  (total  relativistic 
energy),  the  HOMO  (highest  occupied  molecular  orbital 
theory)  energy,  and  the  HOMO-1  (one-level  below  to 
highest  occupied  molecular  orbital)  energy  for  'FeOl+  and 
GdO!+  fragments  in  order  to  evaluate  the  quality  of 
contracted  basis  sets  in  the  representation  of  molecular 
environments.  The  quality  of  these  calculations  was 
compared  to  that  obtained  with  the  basis  sets  (24s  13p), 
(29sl7pl2d),  and  (32s22pl6dl0f). 

TABLE  1  shows  the  TRE,  the  HOMO,  and  the 
HOMO-1  energies  for  the  fragments.  According  to  this 
table,  it  can  be  noted  that  the  TRE,  the  HOMO  and 
HOMO-1  energies  obtained  with  the  contracted  basis  sets 
are  close  to  those  obtained  with  the  uncontracted  basis  sets. 
The  differences  in  the  TRE  are  0.129  and  1.09  hartree  for 
1Fe01+  and  *Gd01+,  respectively.  The  HOMO  energy 
shows  a  difference  of  5. 0x1 0'3  hatree  for  the  two 
fragments.  While  for  the  HOMO-1  energy  the  differences 
are  2.0xl0'3  and  1.0xl0'2hartree. 


Table  1.  TER,  HOMO  energy  and  HOMO-1  energy  (hartree)  for  the  1 FeO I+  and  1  GdO1  +molecuIar  fragments 


Molecular 

fragment 

Basis  sets 

TRE 

HOMO  energy 

HOMO-1  energy 

1Fe01+ 

[12s6p5d]V[4s2p]2 

-1336.314429 

-0.281 

-0.663 

(29sl7pl2d)3/(24sl3p)4 

-1336.443455 

-0.286 

-0.665 

3Gd01+ 

[  1 9s  1 2p8d4f]  5/[4s2p] 2 

-10892.24445 

-0.219 

-0.569 

(32s22pl6dl0f)6/(24sl3p)4 

-10893.33562 

-0.224 

-0.579 

1,3Contracted  and  uncontracted  basis  sets  for  Fe  (5D)  atom;  2,4Contracted  and  uncontracted  basis  sets  for  O  (3P)  atom. 
5,6Contractd  and  uncontracted  basis  sets  for  Gd  (7F)  atom. 


The  properties  of  polyatomic  systems  are  best  basis  sets  used  in  calculations  with  these  systems.  As 

described  when  polarization  functions  are  included  in  the  forGdFeOs,  the  polarization  function  was  considered  in  the 
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[4s2p]  basis  set  for  O  (3P)  atom.  It  was  extracted  from  de 
primitive  basis  set.  The  adequate  exponent  was  chosen 
through  successive  calculations  for  GdFeC>3  by  using 
different  primitive  functions,  taking  into  account  the 
minimum  energy  criterion  at  DKH  second-order  scalar 
relativistic  level.  The  polarization  function  obtained  by  this 
strategy  isad=  0.396928. 

In  order  to  describe  the  configuration  of  a  metal  in 
a  polyatomic  system  it  is  necessary  to  include  diffuse 
functions  in  the  basis  set  for  the  metal.  The  configurations 
of  the  metals  in  GdFeCb  were  adequately  described  by 
adding  a  function  by  symmetry  to  the  basis  set  of  each 
metal  atom.  The  diffuse  functions  obtained  through  the 
total  energy  optimization  of  the  ground  state  anions  Fe1" 
and  Gd1"  by  GCHF  method.  For  the  12s6p5d  and 
19sl2p8d4f  basis  sets,  the  diffuse  functions  are:as  = 
0.0145512;  ap  =  0.115324;  ad  =  0.0538269andas  = 
0.00873780;ap  =  0.0617874;ad  =  0.262367;  af  =  0.118836, 
respectively. 

In  this  work,  atomic  calculations  (contraction  of 
the  basis  sets  and  choice  of  diffuse  functions)  were 
performed  with  the  ATOMSCF  program  [21],  while 
calculations  with  molecular  systems  were  carried  out  with 
the  Gaussian  program  [22], 


III.  RESULTS  AND  DISCUSSION 

Before  starting  the  presentation  and  discussion  of 
the  results  obtained  in  the  investigation  of  piezoelectricity 
inGdFeOa,  objective  of  this  work,  it  is  important  to  note 
some  considerations  about  the  fragment  model  that  was 
used  in  the  representation  of  the  crystalline  system  under 
study. 

Fig.  1  shows  the  fragments  we  have  used  as  a 
model  to  simulate  the  conditions  necessary  to  the  existence 
of  piezoelectricity  in  ABO3  perovskites  [11,  12,  23,  24]. In 
Fig.  1  (A);  (a)  represents  the  [GdFe03]2  fragment  having 
the  Fe  atoms  fixed  in  the  space;  ( b )  represents  the 
[GdFeOsJifragment  in  which  the  Fe  atoms  are  being  moved 
+0.005  A  in  the  symmetry  X  axis,  while  the  Gd  and  O 
atoms  are  maintained  fixed;  (c)  represents  the  [GdFe03]2 
fragment  in  which  the  Fe  atoms  are  moved  -0.005  A  in  the 
symmetry  X  axis  and  Gd  and  O  atoms  are  maintained 
fixed.  Fig.  1  (B)  represents  the  [GdFe03]2  fragment  having 
the  bond  lengths  r(FeiCh),  KFeiCT),  r(Fe204)  shortened  in 
0.005  A. The  same  fragmentmodel  was  used  in  the  study 
developed  with  SmTi03  [11],  YFe03  [12],BaTi03  [23]  and 
also  in  investigations  with  LaFe03[24].  In  this  report,  the 
model  represents  the  crystalline  3D  periodic  GdFeOs 
system.  The  Fe  atom  is  located  in  the  center  of  the 


octahedron,  being  wrapped  up  by  six  O  atoms  arranged  in 
the  reticular  plane  (200)  and  two  Gd  atoms  arranged  in  the 
reticular  plane  (100). 


(A) 


(B) 

Fig.  1.  The  octahedral  [GdFeO 3] 2  fragment.  (A):  (a) 
represents  the[GdFeO 3)2  fragment  having  the  Fe  atoms 
fixed  in  the  space;  (b)  represents  the  [GdFeO 3] 2  fragment 
in  which  the  Fe  atoms  are  being  moved  +0.005  A  in  the 
symmetry  X  axis,  while  the  Gd  and  O  atoms  are 
maintained  fixed;  (c)  represents  the  [GdFeO  3]  2  fragment  in 
which  the  Fe  atoms  are  moved  -0.005  A  symmetry  X  axis 
and  Gd  and  O  atoms  are  maintained  fixed;  (B):  represents 
the[GdFe03]2  fragment  having  the  bond  lengths  fiFejOs), 
fiFeiOfi,  and  r(Fe20fi  shortened  from  0.005  A 


It  is  important  to  note  that  in  order  to  study  the 
crystalline  3D  periodic  GeFeOssystem  it  is  necessary  to 
choose  a  fragment,  a  molecular  model,  capable  of 
adequately  represent  the  physical  property  of  the 
crystalline  system  as  whole. 
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We  still  would  like  to  point  out  three  very 
important  strategic  aspects  in  our  theoretical  approach  in 
the  study  of  gadolinium  ferrite  piezoelectricity,  namely:  (i) 
Firstly,  we  consider  that  the  piezoelectric  properties  in 
GdFeOjresult  from  uncentersymmetric  characteristics 
presented  by  central  ion  and  the  probable  polarization  of 
the  crystal  when  submitted  to  mechanical  stress. (ii) 
Secondly,  the  geometry  optimization  of  the  [GdFe03]2 
fragment  in  the  Cs  symmetry  and  electronic  state  'A  was 
carried  out.  (iii)  Finally,  single-point  calculations  were 
performed  with  the  optimized  geometry,  according  to  the 
descriptions  shown  in  Fig.l,  and  their  results  were 
analyzed  from  the  point  of  view  of  strategy  (i). 


perovskite  GdFe03.Initially,  we  would  like  to  highlight  the 
comparison  of  calculated  and  experimental  bond  length  (r) 
values  [25]  and  then  examine  the  quality  of  the  contracted 
basis  set  supplemented  with  polarization  function  and 
diffuse  functions,  [4s2pld]/[13s7p6d]/[20sl3p9d5f],in 
describing  the  geometry  of  the  system  under  investigation. 
In  TABLE  2,  we  show  and  compare  results  for  r(1Fe01+) 
andr(1GdO!+).  In  this  table,  the  theoretical  results 
reproduce  very  well  the  experimental  values  with 
differences  (A r)  of  4.32  x  10'2Aand  8.00  x  10"2  A, 
respectively,  showing  the  good  performance  of  the 
molecular  [4s2pld]/[13s7p6d]/[20sl3p9d5f]  basis  set  in 
the  description  of  the  geometry  of  the  [GdFe03]2fragment 
studied. 


We  will  now  return  to  the  presentation  and 
discussion  of  the  results  obtained  in  this  study  with  the 

Table  2.  Theoretical  and  experimental  bond  length  (A)  for  lFe01+  and  1 GdO,+  molecularfragments 


Bond  length 

Molecular  fragment 

Theortical 

Experimental  [25] 

Ar 

r('Fe01+) 

1.9063 

1.9495 

4.32  x  lO'2 

r(1Gd1+) 

2.1852 

2.2652 

8.00  x  10'2 

A r  =  ITheoretical  -  Experimental! 


TABLE  3  presents  the  dipole  moment  and  the  TER  energy 
of  the  [GdFe03]2  fragment.  It  can  be  seen  from  this  table 
that  the  result  with  Fe3+central  ion  at  ( b )  or  (c)  is  0.014  and 
0.35  hartree,  respectively,  more  stable  than  at  ( a ),  pointing 
out  that  the  Fe3+  central  ion  is  not  centersymmetric.  This 
evidence  is  confirmed  by  the  decrease  in  dipole  moment  at 
positions  (b)  and  (c)  when  compared  to  the  Fe3+central  ion 
position  at  (a). In  addition,  the  table  also  shows  that  the 

.Table  3.  Dipole  moment  (D)  e  TER 


decrease  in  r(FeiCL),  AFiOr),  and  r(Fe20,t)bond  lengths, 
the  calculations  at  the  atomic  position  (</), provoke  a 
migration  of  atomic  charges  among  the  atoms  of  the 
[GdFe03]2  fragment  and,  consequently,  a  decrease  of  the 
dipole  moment.  These  results  suggest  that  there  is  a 
probable  polarization  of  the  GdFe03  crystal  when 
submitted  to  mechanical  stress. 

rtree)  of  the  [GdFeOj]  2  fragment 


Fe  atom 
position 

Fx 

Py 

pz 

F 

TRE 

(a) 

-0.25 

4.19 

2.08 

7.61 

-24609.25295 

(b) 

-0.613 

-1.86 

1.87 

7.46 

-24609.26693 

(c) 

-2.04 

0.718 

0.0311 

2.16 

-24609.60649 

(d) 

2.99 

3.51 

1.14 

4.75 

-24609.66129 

TABLE  4  shows  the  total  atomic  charges  of  the  [GdFe03]2 
fragment  when  Fe3+  central  ion  is  at  positions  (a),  (b).  and 
(c),  and  for  the  shortened  rlFeiOs),  r(Fei04),  and  r(Fe204) 
bond  lengths,  position  ( d )  .  In  this  table,  when  the  Fe3+  ion 
central  is  moved  to  position  ( b ),  atomic  charges  migrate 
from  the  Fe2,  Gdi,  and  Gd2  atoms  to  the  O2,  O3, 04,  O5,  and 
06  atoms.  On  the  other  hand,  when  the  Fe3+  central  ion  is 
moved  to  position  (c),  a  migration  of  atomic  charges  from 


the  Fe2,  Gdi,  and  Gd2  atoms  to  the02,03,  O5,  and  06  atoms 
also  occurs.  Shortened  AFeiO-s),  AFeiOr),  and 
r(Fe204)bond  lengths,  position  (d)  provoke  a  charge 
migration  from  the  Fe2,  Gdi,  and  Gd2  atoms  to  the  O2,  O5, 
and  06  atoms.  We  can  notice  that,  according  to  TABLE  4, 
the  oscillation  of  the  Fe3+  central  ion  and  the  shortening  of 
the  bond  lengths  caused  by  mechanical  stress  provokes  a 
rearrangement  of  the  total  charges  on  atoms  compared  with 
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the  position  (a)  of  the  [GdFeCblifragment.  The 
rearrangement  ofthe  charges  at  (d)  as  well  as  the  decrease 
of  the  dipole  moment  can  lead  us  to  suppose  that  the 
decrease  of  r(FeiOi),  HFeiOr),  and  rO'ejOr)  bond  lengths 
provokes  a  polarization  of  the  [GdFeOs^  fragment.  We  can 
also  report  that  the  electron  density  does  not  change 
significantly  at  atomic  positions  ( b )  and  (c)  -  the  Fe3+ 
central  ion  is  deviated  from  the  main  axis  of  X  symmetry  - 
when  compared  with  the  atomic  position  (a).  This  is  an 
indication  that  the  nature  of  Fe-O  and  Gd-0  chemical 
bonds  does  not  change  when  the  Fe3+central  ion  deviates 
from  the  X  symmetry  of  the  [GdFe03]2  fragment.  At  (d) 
position,  the  decrease  of  the  r(FeiOi),  rlFeiOr),  and 
r(Fe204)bond  lengths  provokes  the  deviation  of  electronic 
density  in  the  [GdFe03]2  fragment,  indicating  the  nature  of 
the  Fe-O  and  Gd-O  chemical  bonds  was  altered. 


Table  4.  Total  atomic  charges  from  [GdFeO  3]  2  fragment 


Atom 

Atom  position 

(a) 

(b) 

(c) 

(d) 

Flel 

1.60 

1.57 

1.49 

1.48 

Fe2 

0.339 

0.551 

0.571 

0.584 

Gdi 

1.32 

1.68 

1.41 

1.42 

Gd2 

1.21 

1.38 

1.39 

1.43 

Oi 

-1.06 

-0.652 

-0.974 

-0.979 

O2 

-1.32 

-1.34 

-1.49 

-1.49 

O3 

-0.247 

-0.267 

-0.284 

-0.168 

04 

-0.416 

-0.779 

-0.122 

-0.359 

O5 

-0.968 

-0.973 

-0.984 

-0.97 

o6 

-0.864 

-0.870 

-0.953 

-0.949 

Thus,  we  can  notice  that  electrostatic  interactions 
play  an  important  role  in  the  constitution  of  the  electronic 
structure  of  the  [GdFe03]2  fragment.  This  is  consistent  for 
two  reasons,  namely:  (1)  due  to  the  repulsive  effect  of  d 
and  f  electrons  in  both  high-spin  and  low-spin  octahedral 
species  of  ML  complexes  (M=metal  and  L=ligant),  all  d 
and  f  electrons  density  will  repel  the  bonding  electrons 
density  [26];  (2)  in  lanthanide  complexes,  chemical  bonds 
of  ionic  nature  predominate  [27], 

Analysis  of  the  oscillation  of  the  Fe3+  central  ion 
between  the  positions  ( b )  and  (c),  and  the  shortening  of  the 
r(Fei03),  r(Fel04),  and  r(Fe204)  bond  lengths  in  the 
[GdFe03]2  fragment,  position  (d),  showed  that  the 
occurrence  of  the  uncentersymmetrical-centered  cation 
(Fe3+)  and  the  polarization  of  the 
[GdFe03]2fragmentproduce  an  electric  current,  leading  us 


to  infer  in  this  work  evidence  that  GdFe03  may  present 
piezoelectricity  caused  by  electrostatic  interactions  of  its 
atoms. 

IV.  CONCLUDING  REMARKS 

The  use  of  DKH  second-order  scalar  relativistic 
approach  together  with  qualified  basis  sets  to  represent  the 
GdFe03  polyatomic  system  allowed  the  simulation  of  the 
conditions  necessary  for  the  existence  of  piezoelectricity  in 
this  material.  The  calculated  properties  allow  us  to  infer 
that: 

1.  The  Fe3+  central  ion  in  the  [GdFe03]2  fragment  has 
uncentersymmetrical  characteristics.  In  addition,  when 
submitted  to  mechanical  stress,  the  [GdFe03]2  fragment 
presents  polarization,  which  leads  us  to  suggest  the 
existence  of  piezoelectricity  in  GdFeOs,  caused  by 
electrostatic  interactions. 

2.  The  results  of  the  computational  investigation  presented 
in  this  work,  punctuating  the  possible  piezoelectricity  in 
GdFeOs  through  the  use  of  the  fragment  model  used  to 
investigate  this  property  in  BaTi03,  LaTi03,  and  SmTi03 
corroborate  the  viability  of  the  model  and  methodology  to 
investigate  possible  piezoelectricity  in  other  ABO3 
perovskites. 
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